On two descriptions of disordered phase of 2D quantum antiferromagnets: monopole 

plasma and Kalb-Ramond fields 
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We have studied a system of light bosonic spinons interacting with compact gauge fields. By the 
generalization of the works by Polyakov on compact gauge fields the system is mapped to monopole 
plasma model and a model of open surfaces coupled to antisymmetric Kalb-Ramond gauge fields. 
The monopole correlation function in the presence of light spinon is computed based on these two 
models. 



The quantum critical phenomena and the associated 
quantum phase transitions are one of the most outstand- 
ing problems of current physics Q, 0] • These problems 
have been addressed mostly in the context of quantum 
magnetism and heavy fermion systems Q. The stan- 
dard framework for the understanding of general crit- 
ical phenomena is the Landau-Ginzburg-Wilson(LGW) 
theorydH. In LGW theory, critical phenomena are de- 
scribed in terms of the order parameters characterizing 
either side of transition. 

Recently a new framework of quantum criticality be- 
yond LGW theory was proposed 111 ■ This framework 
predicts degrees of freedom emerging at critical point 
which can not be expressed by order parameters of phases 
surrounding the critical point. A very rough picture is 
that the emerging degrees of freedom are confined away 
from critical point due to strong singular ( or topolog- 
ically nontrivial) gauge field fluctuations, while at the 
critical point the singular fluctuations become irrelevant 
(in the sense of renormalization group) liberating the con- 
fined degrees of freedom. As quarks cannot be described 
by mesons or baryons at low energy, the emerging decon- 
fined degrees of freedom cannot be described by order 
parameters of either phase of critical point. 

The key step of the above scenario is to show the ir- 
relevancy of singular gauge field fluctuation at critical 
point. This point has also been investigated intensively 
in the context of massless (2+l)-dimensional quantum 
electrodynamics along with some controversies 0, H, ■ 
In (2+1)- dimensional system the singular gauge field 
fluctuations can be identified with monopole field config- 
urations [lfl Hll | . The monopoles can exist only if the 
associated U(l) gauge field is an angular variable ( in 
other words, has its value on cirle S* 1 ), and such gauge 
field is called compact. Ordinary U(l) gauge field has 
its value on R 1 and is called non-compact. Polyakov 
showed that the static and heavy charges interacting via 
compact gauge fields have linearly rising interaction po- 
tential energy (V(R) oc R), thus they are in confinement 
|TH. T ne potential energy V(R) is the energy of bundle 
of electric flux line between two charges, so it is natu- 
rally proportional to the length of flux line. The result 
by Polyakov [Jjj cannot be directly applied to the prob- 



lem of quantum criticality since the matter fields with 
gauge charges are dynamic and light in the vicinity of 
quantum critical point. 

In this paper we study the descriptions of the disor- 
dered phase two-dimensional quauntum antiferromagnets 
where the interaction between the light bosonic matter 
degrees of freeom and the compact gauge field (especially 
monopole configurations) can be expressed in a very illu- 
minating form. This formulation enables us to compute 
the correlation functions of monopole operators in the 
presence of light spinon. 

We start with the 0(3) nonlinear sigma model formu- 
lation of spin-S* two-dimensional quantum antiferromag- 
nets 2j. 

S n = Sq + Sb 

So = Yg J dT J Sr [ (5rA)2 + 

Sb = e ? I dr I dw n • <9 r n x <9 u n, (1) 

f Jo Jo 

where n = n(r, r) is a three component Neel order pa- 
rameter and ep = (— i) r =+ r H . Sb is the Berry phase term. 
f is to be understood as both continuum and lattice co- 
ordinates. The large N expansion study of So shows that 
the Neel ground state is realized for g 2 < g 2 r , while disor- 
dered massive uniform quantum paramagnets obtains for 
the strong coupling regime g 2 > g 2 r 0] . Haldane showed 
that the Berry phase term is non-vanishing only for the 
monopole configuration of Neel order parameter |l2l | . 

The monopoles can be described by U[l) gauge field 
in the CP 1 representation of Eq. 
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h = Z^crZ, J2\ z a\ 2 = 1 

a=l 



F, 



n V — - n ■ 9 M n x <9„n. 



(2) 



Z is called bosonic spinon in the context of quantum an- 
tiferromagnetism. The compact nature of gauge field 
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is manifest since it is the fiber of Hopf bundle S 3 (Z) — 
S 2 (n) . The monopole of strength p is characterized by 



of an operator 



Trln 



±- / F, F = ^-F^dx^ A dx u , (3) = i [°° ^L e -rm%/2 f 

2lr J* 2 2 J e t A 



(-ia - a) 2 + ?n 



dr. 



-Tre"2 



[(—id—A)+m^] 



where E is a closed surface enclosing the monopole. Read 
and Sachdev showed that the Berry phase term can be 
reduced to [3 



S B = -i(2S)~£ CsPs, 



(4) 



where s denotes the site of dual lattice and p s is the 
strength of monopole located at s. Cs = 0, 1, 2, 3 is an in- 
teger valued field designating sublattice structure of dual 
lattice. We note that for p s = 4, e~ SB = 1, so that the 
Berry phase term disappears. 

The CP 1 model can be naturally generalized to CP 
model and it possesses a critical point. Relaxing the con- 
straint Y2a=i \Za\ 2 — 1 by mass term mz, we obtain a 
model to be investigated in this paper. 



S = Sz + Sa + Sb 
Sz = J d 3 x - iA^)Z\ 2 + m 2 z \Z\ 2 ) 



(5) 



where Sb is given by Eq. Q). e 2 (which has dimension 
of mass) will be treated as phenomenological parameter 
independent from mz- This model is a slight modifi- 
cation of the one considered by Sachdev, Jalabert, and 
Park [H 03. In their works the singular monopole gauge 
fields are expressed separately in action. We will assume 
that mz is light enough, so that only the monopoles with 
strength p = 4 mod 4 survive in the functional integral. 
If mz is large enough the matter Z fields can be simply 
dropped in low energy sector. The monopole summation 
(see below for details ) can be done with Berry phase 
taken into account, and it leads to cos[(^ s — 'TrS'Cs]. The 
sine-Gordon interaction pins the tp s , thus depending on 
values of S spin-Peierls states with broken lattice trans- 
lational symmetry obtain 

We can formally integrate out the spinon field Z ex- 
actly from Eq. ©. It contributes 



5S z = NTi\n (—idp — A^) 2 



to the effective action of gauge fields. Next we employ 
the well-known path integral representation of logarithm 



x exp 
i 



X(0)=X(r) 



= ^W(C), W(C) 



Jc 



C 



(6) 

C denotes the closed trajectory of the Z-particle. W(C) 
is the Wilson loop [1^ . e is a short time cut-off. 
Then the partition function reads 

Z = J D[A]e- s *- N X'<> w(c) 

= £ u z r - £-£w c o---^(c ro )>-(7) 



to! 

m=0 Ci C m 

The average of Wilson loop is given by 



(W(Ci) ■ ■ ■ W(C m )) = I D[A] e- s ^J e E > §c > A 

' D[A)e- s ^e iJ: ^^ F , (8) 



where T,j is an open surface whose boundary is Cj. The 
gauge field can be decomposed into non-compact and 
monopole part 10]. Accordingly the averages of Wilson 
loops are factorized 

(W(Ci) • • • W(C m )) = {W{C X ) ■ ■ ■ W(C m )) 

x{W{C 1 )---W{C m )) M . (9) 

At low energy the monopole contributions dominate and 
we focus on the {W(C X ) ■ ■ ■ W{C m )) M - 

The partition function of monopole plasma is given by 

nana 



JV 



Za = £ m / II dx * ex P 

JV=0[{g o } ■ J 3=1 



TV 



' 2e 2 ^ |x Q 



QaQb 



The fugacity is 



= A 3 



,-ciA/e 



Xfc 

(10) 
(11) 



where A is a high energy cutoff and c\ is a numerical 
constant of order unity. The summation over monopoles 
can be done with Wilson loops inserted, and it results in 

(W(Ci) ■ ■ ■ W{C m )) M = ^r I D[<p] 



z a 



x exp 



g-j / d 3 x(Wf) 2 + 2z I d 3 xcos[p(<p + n)} 



Z = I D[<p] exp - — / d 3 x(Wip) 2 + 2z I d 3 xo 
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where 



1 (y 3 - x) • da( yj ) 



(13) 



is the source term of multiple Wilson loops, yj is the 
coordinates on the surface Sj. Eq. I|13|) is nothing but 
the solid anlge subtended by surface X observed from 
the point x. If particles carry nontrivial spinor, color, 
or flavor indices, the mutiple Wilson loops contain non- 
commuting matrices, so that it cannot be expressed in a 
simple way as Eq. 1|12|) . f represents the photon in dual 
representation. As can be seen in Z of Eq. (|12|) (p is 
pinned at the bottom of cosine potential and the photon 
becomes massive with exponentially small mass 



A 



A\3/4 



-ciA/2e 



(14) 



The partition function Eq. and the average of Wilson 
loops of the form Eq. I112H with appropriate operator 
insertions constitute a monopole plasma representation 
of light bosonic matter field coupled to compact gauge 
fields. 

Doing the monopole functional integral via steepest 
descent for a single large specified loop we find that 



(W(C)} ~ exp[-CT5] 



e 2 m r 



(15) 



where S is the minimal area bounding the trajectory C 
|llj | . From the viewpoint of Hamiltonian formulation Eq. 
(I15|l implies the linear potential, thus confinement. 

In confining phase there exist electric field flux lines 
connecting two charges, and the time evolution of these 
lines generate open surfaces. Therefore, it is tempting to 
attempt to express the average of Wilson loop Eq. H12J1 in 
terms of these surfaces. This has been done by Polyakov 
(lif. Consider a functional integral 

Z[{W(C)}}= J D[0,B]e- r ^ B ^ 

T[cf>,B,Ec} = J d 3 x 

+ (2z)(l - cos p(f>) 



— B^B^ + — <\>t ( !B 



i 

2tt' 



i l B. 



(16) 



B^y is an antisymmetric Kalb-Ramond gauge field 0] 
which naturally couples to surface just as vector potential 
couples to point particle in a minimal way J dx^A^. B — 
^B^dx^ A dx v is a differential 2-form. £<7 is an open 
surface bounding closed path C. The generalization to 
multiple surfaces is immediate. 



i I B 

'So 
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B. 



With B^y eliminated by solving equation of motion, it 
can be shown that 

W (C)W M (C) <x Z[{W(C)}\, <j>-r, = tp. (17) 



This establishes the equivalence between monopole 
plasma approach and the Kalb-Ramond gaug field ap- 
proach. Now one can integrate over <p first to obtain 
effective action of B^ v , and the result still should be the 
same as the Wilson loop W{C). In three dimension 

dB = Hdx 1 A dx 2 A dx 3 , H = -e^dxB^. 

Eliminating </> by solving equation of motion of 4> we ob- 
tain 



W({Cj}) 



D\B] e- S ^ 



-if* 



S(B) = J d 3 x(-L B ^ 



1 

2np 



H sin 



H 



H 



H 2 



where 



m# = A 



A\i/6 



-ciA/3e 



,(18) 



(19) 



Polyakov has obtained a non-standard string theory by 
summing over branches of sin -1 (lij. The partition 

function Eq. J7J) and the averages of Wilson loop Eq. I|18|) 
constitute a description in terms of Kalb-Ramond anti- 
symmetric gauge field. Picking up a branch near H = 
and expanding in H for large loops the effective action 
S(B) becomes 



i 



(V • Bf 



where B 1 — ^e t: > k Bjk- The propagator is 



(B l (x)B 3 (y))=6 l 



pm 3 H e mp l x y l 



(20) 



(21) 



Now the averages of Wilson loop can be computed from 
Eq. CHI) 



W{{Cj}) = exp 



pm 3 H 



E 



da a ( Xi )da a (yj) 



e -™p ! x i-yj 
4vr|x l -y i | 



(22) 



For a very large loop the main contributions will come 



from the region |x — y| ~ m p 1 , da a 
tation value of single loop is then 



The expec- 



(W(C)} ~ exp - m pCl e 2 / d 2 ^ , (23) 



where 



g a b = d a x ■ d b x (24) 
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is the induced metric. Thus the minimal surface with 
prescribed boundary will give the dominant contribution 
to the Wilson loop (see Eq. dJ). 

Now let us discuss the correlation function of monopole 
operator ( e W« e -w(")) with |z - w| > m" 1 . 

We note that the monopole operator can be identified 
with the valence bond operator 0, [2(| • 



Jip(r) 



(25) 



Since the monopole operator acts trivially on non- 
compact sector the correlation function can be expressed 
by 



/giP¥>(z)g-w(w)\ 



— E 



x (e lpv{z) e 



{-l) m N r< 



m i 



X\(W{C 3 )) M , 



(26) 



where Zm is the partition function in monopole sector 
{-l) m N n 



Zm 



E 



-^■■■Z'cJllWiC^M. (27) 



Let us calculate m = and m = 1 contributions to the 
correlation functions. The spinon is absent in m = sec- 
tor, and the correlation function is essentially determined 
by pinning and massive photons. 



( e w(z)g-w(w)^ m=o 



1 



;/l„|z-w| 



(28) 



In the computation of m = 1 contribution it is essential 
to include the m = 1 part of Zm- 



Z M (m = l) = -NJ2(W(C)) M 

c 



-N 



i 

E 

c 



-N 



exp[-aS(C)] 

3 — T -D/2 -m%T/2-aS(C) 



(29) 



where the second line is from Eq. I|15|l and the summa- 
tion over closed paths can be done by following Chap. 9 
of Ref . |Toj| . In the above mz is to be understood as a 
mass in continuum theory. I? is a space-time dimension. 
For spinless particle the length of typical dynamical path 
r is the order of area S (C) (like Brownian motion) [2l| . 
(For fermions r ~ sj S(C)) 

Next we compute the m — 1 part of the numerator of 
Eq.®- 



1 



- N 
A) 



D[ip] exp 

-2zcos(p(ip + 1])) + ipip(w) — ipip(z) 



(30) 



For a large loop the functional integral can be computed 
in a saddle point approximation. (</> = ip + r/) 



— V 2 (0 — rf) = 2pzsinptf>+ipS(x—z) — ipS(x. — w). (31) 

4-7T 

V 2 r] is the source term from the spinon loop. 

The investigation of Eq. I|31() reveals that when z, w 
are far away from surface bounding boson loop, the cor- 
relation is rather small. 



exp 



-(75(C) 



-m n z— w 



Main contributions come from the cases where the 
monopole operators are located on the surface Ec- ln 
this case the local perturbation of ip field can be relaxed 
by soft transverse modes of surface whoe energy is of the 
order 



E 



1 



(32) 



where R is the size of the closed path ^ij. This is the 
way how matther fields interact with monopoles. Thus 
for the surfaces on which the monopole operators reside 
the correlation becomes 



exp 



tS 



-|z-w|/e 2 fl 2 



Taking into the normalizations Eq. Q29fl we obtain 

c 



/ e w(z)g-w(w) 





"g-|z-w|/S" 




(exp 








. |z-w| _ 





— I ; % ' p -tS-\-w-z\/s 

|w — z, 

1 1 c 



(33) 



where S(C) ~ R 2 is used. Employing the scaling r 
S(C) 



^ e -<rS-|w-z|/S 



C 

x exp 



dr 



- 7 /2 



4' 



m 7 j w — z 
— ut 



mt + a 



-7/4 



+ (34) 



where K\ is the modified Bessel function. From Eq. H34|) 
it is clear the massive boson suppressess the monopole 
fluctuation. The exponent 7 cannot be put to 3 since only 
the subclass of surfaces contribute to Eq. Since the 

addition of one more point determines a plane 7 may be 
approximated to be 1. The monopole correlation func- 
tion is significantly influenced only in the limit where the 
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renormalized spinon mass ro| + 2er becomes exceedingly 
small. In this limit 



const. 



17/4+3/2 ■ 



(35) 



The result Eq. I|35|l clearly shows that light spinons tend 
to prevent the proliferation of monopoles. The contribu- 
tions from m > 2 terms involve various configurations of 
surfaces, such as foldings and intersecting surfaces. This 
problems are left for future investigations. 

In summary we have investigated two descriptions of 
the system of bosonic spinons interacting with compact 
gauge field. These descriptions provide rather clear pic- 
ture on the interaction of light spinons with monopole 
excitation of gauge field. Monopole correlation functions 
is computed up to the first order of spinon loop and it 
turns out that the soft modes of the surface bounding 
spinon trajectory plays a crucial role. 

The author is thankful to Chanju Kim for useful com- 
ments. 
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